Magnetism in Geometrically Frustrated YMnC>3 Under Hydrostatic Pressure Studied 

with Implanted Muons 
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The ferroelectromagnet YM11O3 consists of weakly coupled triangular layers of S = 2 spins. 
Below Tn « 70 K muon-spin relaxation data show two oscillatory relaxing signals due to magnetic 
order, with no purely relaxing signals resolvable (which would require different coexisting spin 
distributions). The transition temperature Tn increases with applied hydrostatic pressure, even 
though the ordered moment decreases. These results suggest that pressure increases both the 
exchange coupling between the layers and the frustration within the layers. 

PACS numbers: 75.50.Ee, 76.75. +i, 75.40.Cx, 75.47.Lx 



The hexagonal i?Mn03 manganites d (R=Eo, Er, 
Tm, Yb, Lu, Y or Sc) are a class of magnetically or- 
dered materials that also possess ferroelectric properties 
and some degree of magnetoelectric coupling, suggesting 
that an understanding of the magnetism may allow ma- 
nipulation of the electric polarization with possible device 
applications 0. Moreover, these compounds are layered 
and exhibit geometric frustration within their layers, of- 
fering the possibility of studying the effect of competing 
interactions in low-dimensional systems Q. 

YMnC>3, an insulator which undergoes a ferroelectric 
transition at Te = 913 K [4j, is the most intensively 
studied of the i?MnC>3 series. The magnetic system is 
based on a frustrated architecture, with Mn 3+ (S = 2) 
ions forming a two-dimensional (2D) corner-sharing tri- 
angular network. Studies of magnetic susceptibility [B} 
confirm the frustrated nature of the system, with a large 
ratio of Weiss temperature (|0| = 705 K) to antiferro- 
magnetic (AFM) ordering temperature (Tn « 70 K). 
Neutron diffraction studies 0, 0] have found that be- 
low Tn, the Mn 3+ spins lie in the ab plane and adopt 
a 120° structure, with a Mn moment of 2.9-3. 1/ib at 
1.7 K. This is below the expected value of 4/iB, due to 
fluctuations associated with either the frustration or the 
low dimensionality. Heat capacity measurements Q 
suggested incomplete ordering of the Mn spins below Tn 
and elastic and inelastic neutron scattering (INS) mea- 
surements @ found strong diffuse scattering persisting 
across T N ; this was taken as evidence for a spin liquid 
phase which coexists with the ordered phase below Tn- 
However, it has been claimed that by taking low energy 
Einstein modes into account, heat capacity data may be 
consistent with conventional AFM ordering [lOj. Other 
INS studies [ll|, [r| have found evidence for coexisting 
three-dimensional (3D) and 2D fluctuations. More recent 




FIG. 1: ZF At+SR spectra measured at T = 6.1 K and 80.4 K 
at p = 1 bar, with oscillations clearly observable below Tn. 



neutron studies 13, UH show that the ordered moment 
decreases with increasing hydrostatic pressure. This ob- 
servation was explained in terms of a pressure-induced 
change in volume fraction of ordered and spin liquid com- 
ponents of the material jl 31 ] - 

In this paper we present the results of muon-spin re- 
laxation [l5j (/U + SR) measurements made on YMnC>3 at 
ambient pressure and as a function of hydrostatic pres- 
sure p up to p = 13.7 kbar. Muons are a sensitive local 
probe of the spin distributions in a magnetic material 
and have proven particularly useful in probing frustra- 
tion related effects 16j; however, measurements made 
under pressure have been less common. Strikingly, we 
show that Tn increases with increasing p, even though 
the magnetic moment decreases. This provides, as we 
shall show, an insight into the role of the finely balanced 
interactions in this system. 

Zero- field (ZF) /i + SR measurements [TtJ were made on 
a single crystal sample of YMnC>3 using the GPS instru- 
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ment at the Swiss Muon Source (S/iS) and on the MuSR 
instrument at the ISIS facility. The muon spin was di- 
rected along the crystallographic c-axis of the sample. 
Transverse field (TF) /x + SR measurements were made on 
a powder sample of YMnC>3 under hydrostatic pressure 
using high energy incident muons on the /iEl decay muon 
beam line at S/xS. For these measurements, the sample 
was packed into a cylinder approximately 7 mm in diam- 
eter and 18 mm in length which was mounted in a Cu- 
Be25 piston cylinder pressure-cell with Daphne oil used 
as the pressure medium. The pressure was measured by 
monitoring the superconducting transition temperature 
of an indium wire located in the sample space [lj| . 

Example ZF /i + SR spectra at ambient pressure (mea- 
sured at S//S) are shown in Fig. Q] Below Tn we observe 
oscillations in the time dependence of the muon polar- 
ization (the "asymmetry" [l5j) characteristic of a quasi- 
static local magnetic field at the muon stopping site. This 
local field causes a coherent precession of the spins of 
those muons with a component of their spin polarization 
perpendicular to this local field. The frequency of the 
oscillations is given by v$ — "f f j i Bi/2n, where 7 M is the 
muon gyromagnetic ratio (= 2ir x 135.5 MHz T _1 ) and 
Bi is the average magnitude of the local magnetic field at 
the zth muon site. Any fluctuation in magnitude of these 
fields will result in a relaxation of the signal, described 
by relaxation rates A^. Two separate frequencies were 
identified in the low temperature spectra, corresponding 
to two magnetically inequivalent muon stopping sites in 
the material. The larger frequency v\ is found to have 
a small relaxation rate Ai while the smaller frequency 
i>2 is associated with a relaxation rate A2 which is an 
order of magnitude larger. The spectra were found to 
be well fitted using only oscillatory components. In the 
ordered phase of YM11O3 the Mn 3+ moments adopt the 
120° structure, where the ordered Mn moments lie within 
the a-b plane [6| . Although the initial muon polarization 
is directed parallel to the c-direction, the fact that the 
muon couples to dipole fields means that, in addition to 
the magnetic field components directed perpendicular to 
the muon spin, there may also exist components parallel 
to the muon spin. In our measurements these compo- 
nents only give rise to a constant background offset. 

To follow the temperature evolution of the observed 
features the S^iS spectra below Tn were fitted to the func- 
tional form 

2 

A(t) = A hg + ^e A,t cos(2tu^ + (pi), (1) 

i=l 

where A\, g represents a constant background contribu- 
tion, including the signal from those muons that stop in 
the silver sample holder or cryostat tails. Nonzero phases 
(pi were required to fit the observed oscillations because 
of the difficulty in resolving features at early times in 
the spectra due, at least in part, to the fast initial de- 
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FIG. 2: (a) Relaxing asymmetry (normalised) A rB \ (measured 
at ISIS) showing a step- like change at Tn (see text), (b) 
and (c) Fits of the S/iS spectra at p = 1 bar to Eq[T] (b) 
Temperature dependence of precession frequencies u\ (open 
circles) and (closed circles). The solid line is a fit to Vi{T) = 
fi(0)(l — (T/Tn)")' 3 (see main text), (c) Relaxation rates Ai 
(open circles) and A2 (closed circles). 



polarization feature. To give the best fit across the en- 
tire temperature range we fixed A\ = 4.0%, (pi — —34° 
A 2 = 8.8% and 2 = -19°. 

We are unable to resolve a relaxing signal due to any 
coexisting spin distribution such as the spin liquid phase 
suggested to persist into the ordered phase [9( . Our muon 
data are straightforwardly accounted for by a model as- 
suming conventional long range magnetic order through- 
out the bulk of the sample. The full muon asymmetry is 
observed above Tn and relaxes exponentially with a sin- 
gle relaxation rate, as expected for a conventional param- 
agnetic state. Below Tn we observe an oscillatory signal 
with only weak relaxation, as expected for a well-defined 
magnetically ordered state. Thus we find no evidence for 
any static magnetic inhomogeneities in our sample, nor 
any evidence for coexisting ordered and disordered vol- 
ume fractions. Furthermore, simulations of the internal 
field distribution expected from "droplets" of spin liquid 
dispersed in an antiferromagnetic medium predict a size- 
able slowly relaxing fraction, inconsistent with our data. 
This effectively rules out a model in which long range 
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order and a disordered spin liquid somehow coexist in 
different regions of the sample. Thus the diffuse scatter- 
ing observed in neutron experiments Q (which has been 
proposed to arise from a spin liquid state) results not 
from phase separated regions but more likely from high 
frequency fluctuations which are motionally narrowed on 
the muon timescale and which do not affect the mag- 
netic ground state of the system. Our data do show a 
small missing fraction of relaxing asymmetry (approxi- 
mately 10% of the total signal) corresponding to muons 
which are depolarized within 1 ns of implantation (not 
resolvable in our measurements) but this feature only ap- 
pears below Tn, unlike the diffuse scattering which per- 
sists across the transition [9] . 

Complementary data were measured over a longer time 
window at the ISIS facility, where the limited time reso- 
lution does not allow us to resolve oscillations. Instead, 
we see a sharp decrease in the relaxing amplitude A le i as 
the material is cooled through Tn. This is because the 
local field in the ordered state will strongly depolarize 
the muon-spin if local field components are perpendicu- 
lar to the initial muon-spin polarization (or have no effect 
on it if local field components are parallel to the initial 
muon-spin polarization) , removing the relaxing asymme- 
try from the spectrum. These data (Fig. [2(a)) show that 
the relaxing amplitude does not vary with temperature 
apart from at Tn, showing that there is no temperature 
variation in the volume fractions due to the ordered mag- 
netic state (for T < Tn) or paramagnetic state (T > Tn). 
Moreover, only a weak, temperature independent relax- 
ation (« 0.015 MHz) is observed in the ISIS data below 
Tn, which is well within the ordinarily expected back- 
ground contribution, confirming that no purely relaxing 
component is needed in Eq. 

Fig Ob) shows that vi and v% can be fitted by 
Vl (T) = v l {{)){\ - (T/T N ) Q )' 3 from which we estimate 
T N = 74.7(3) K, a « 2.5 and (3 = 0.35(3), consistent with 
3D Heisenberg or 3D XY behaviour, as found in specific 
heat studies [10(. Our determination of Vi(0) allows us 
to attempt to identify the muon sites in YM11O3 . Dipole 
fields were calculated in a sphere containing ~ 10 5 Mn 
ions with moments of 2.9/xb arranged in the 120° struc- 
ture. The positive muon's position is usually in the vicin- 
ity of electronegative O 2- ions [HI. Candidate muon 
sites giving rise to the higher oscillation frequency v\ are 
found to be separated from an 0(4) oxygen by 1 A along 
the c-direction. This gives sites at coordinates (1/3, 2/3, 
z) and (2/3, 1/3, z), where z « 0.09, 0.42, 0.59 and 
0.92. Several candidate sites for the lower frequency v-i 
are found close to the planes of triangularly arranged 
0(1) and 0(2) atoms. One possibility is for the sites to 
lie between oxygens, with the sites forming triangles cen- 
tred again on the (1/3, 2/3, z) and (1/3, 2/3, z) positions, 
where now z f» 0.18 and ~ 0.32. 

The division of the muon sites into a set lying close to 
the Mn planes and a set between these planes may ex- 
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FIG. 3: (Color online.) Results of TF /i + SR measurements 
made under hydrostatic pressure, (a) Example spectra mea- 
sured at p = 11.7 kbar in a TF of 5 mT at T = 71 K (red 
circles) and T = 80 K (black circles), (b) Temperature depen- 
dence of the normalised amplitude, A le i, for several applied 
pressures, (c) The transition temperature Tn is seen to in- 
crease linearly with increasing pressure. 



plain the difference in the observed relaxation rates. The 
relaxation rates are expected to vary as A ~ A 2 r, where 
A is the second moment of the local magnetic field distri- 
bution and r is its fluctuation time. Sites giving rise to 
the frequency V\ lie close to the Mn planes in well defined 
positions such that A is small. We also expect muons at 
these sites to be sensitive to both in-plane (2D) and out- 
of-plane (3D) magnetic fluctuations. Relaxation rate Ai 
increases as Tn is approached from below (Fig. O^c)) be- 
cause r increases due to the onset of critical fluctuations 
close to the phase transition. Sites associated with fre- 
quency vi lie between the Mn planes in several positions 
where there is some variation of the dipole fields around 
40 MHz and consequently a large value of A. These will 
be less sensitive to 2D fluctuations than those lying close 
to the planes, reducing the influence of any variation in r. 
The temperature evolution of A2 (Fig. [2(c)) is now dom- 
inated by the magnitude of A, which scales with the size 
of the local field. Relaxation rate A2 therefore decreases 
as the magnetic transition is approached from below. 



4 



In contrast to the ZF data measured at ambient pres- 
sure at S/zS, it is not possible to resolve precession fre- 
quencies in S/LiS ZF data measured under hydrostatic 
pressure. This is because the signal from the sample 
amounts to only 20% of the total measured spectrum 
from the pressure cell and is strongly depolarized. To 
observe the magnetic transition, it was necessary to fol- 
low the amplitude of the muon precession in a transverse 
field of B t = 5 mT (Fig.[3ja)) and fit the data to the func- 
tional form A(t) = [A + (A^ - A )A Icl ]e- At cos^TJ^), 
where Aq and Aoo are found to be independent of pres- 
sure, A is a relaxation rate and A rc \ is the normalized 
amplitude. The temperature dependence of A vc \ charac- 
terizes the magnetic transition and is shown in Fig. E^b) 
for several applied pressures. The pressure indepen- 
dence of Aq and A^ are inconsistent with an earlier 
speculation that there is a change in volume frac- 
tion between spatially separated spin-liquid and AFM 
states when the pressure is varied. To extract the tran- 
sition temperature we use the phcnomenological form 
A lc \ = [1 + exp(— b{T — Tn})] -1 where b is a parame- 
ter describing the "width" of the transition. At all pres- 
sures the parameter b remained practically unchanged 
at b w 1.2 K _1 , while Tn shows a linear increase with 
pressure as shown in Fig. Etc). A straight line fit yields 
Tn = 74.0 + 0.29p, where p is the pressure in kbar. Our 
value for dT^/dp = 0.29 Kkbar -1 is in excess of that 
predicted by the "10/3" law [20], dT N /dp = (10/3)T N /B 
(where B is the bulk modulus), which holds for many 
oxides and garnets. Using B = 1.65 Mbar [2lJ yields 
an estimate dT^/dp = 0.14 Kkbar -1 , which is approx- 
imately half the measured value. This strong pressure 
dependence of Tn demonstrates the sensitivity of Tn to 
pressure-induced small changes in intralayer J and inter- 
layer coupling J' (see e.g. Ref. liih. 

The increase of Tn with pressure in YMnC>3 is surpris- 
ing given that previous studies [3] have shown that 
the magnitude of the ordered moment decreases with ap- 
plied pressure, implying an increase in the spin fluctua- 
tions that reduce the value of the magnetic moment found 
at ambient pressure. It is plausible that the application of 
hydrostatic pressure to the polycrystalline material has 
two effects. The first is on the structure of the trian- 
gular MnO planes as shown by neutron measurements 
|13j |. At ambient pressures the Mn-0(3)-Mn and Mn- 
0(4)-Mn bond angles and length differ slightly. This re- 
lieves the magnetic frustration to an extent. The neutron 
diffraction measurements show that upon the application 
of pressure the Mn-0(3)-Mn and Mn-0(4)-Mn bond an- 
gles and lengths approach each other. These effects act 
to make the MnO planes more perfect realisations of a 
triangular lattice, causing the exchange coupling along 
the triangular bonds to become more similar as pressure 
is increased. This increased frustration has been pro- 
posed as an explanation for the reduction in magnetic 



moment [HQ. The second effect is an increase in both 
J' and J, which have an exponentially sensitive depen- 
dence on bond distance, making the relative effect on J' 
larger. INS measurements at ambient pressure show that 
J'/J ~ 10" 2 fnl, which is consistent with T N /J ~ 0.5 
(J w 3meV [9j,[llJ]). Pressure will therefore increase J'/J 
and hence Tn- It is clear that in order for this dual ef- 
fect to occur there exists a delicate balance of competing 
interactions in this system. 
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